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The spectral function of the closed-shell Neon atom is computed by expanding the electron self-energy 
through a set of Faddeev equations. This method describes the coupling of single-particle degrees of free- 
dom with correlated two-electron, two-hole, and electron-hole pairs. The excitation spectra are obtained using 
the Random Phase Approximation, rather than the Tamm-Dancoff framework employed in the third-order al- 
gebraic diagrammatic contruction [ADC(3)] method. The difference between these two approaches is studied, 
as well as the interplay between ladder and ring diagrams in the self-energy. Satisfactory results are obtained 
for the ionization energies as well as the energy of the ground state with the Faddeev-RPA scheme that is also 
appropriate for the high-density electron gas. 



PACS numbers: 31.10.+z,31.15.Ar 

I. INTRODUCTION 

Ab initio treatments of electronic systems become unwork- 
able for sufficiently complex systems. On the other hand, 
the Kohn-Sham formulation QJJ] of density functional the- 
ory (DFT) [2] incorporates many -body correlations (beyond 
Hartree-Fock), while only single-particle (sp) equations must 
be solved. Due to this simplicity DFT is the only feasible 
approach in some modern applications of electronic structure 
theory. There is therefore a continuing interest both in devel- 
oping new and more accurate functionals and in studying con- 
ceptual improvements and extensions to the DFT framework. 
In particular it is found that DFT can handle short-range inter- 
electronic correlations quite well, while there is room for im- 
provements in the description of long-range (van der Waals) 
forces and dissociation processes. 

Microscopic theories offer some guidance in the devel- 
opment of extensions to DFT. Orbital dependent function- 
als can be constructed using many-body perturbation theory 
(MBPT) yl 0]. More recently, the development of general 
ab initio DFT fH 0] addressed the lack of a systematic im- 
provement in DFT methods. In this approach one considers 
an expansion of the exact ground-state wave function (e.g., 
MBPT or coupled cluster) from a chosen reference determi- 
nant. Requiring that the correction to the density vanishes at a 
certain level of perturbation theory allows one to construct the 
corresponding approximation to the Kohn-Sham potential. 

A different route has been proposed in Ref. [7] by devel- 
oping a quasi-particle (QP)-DFT formalism. In the QP-DFT 
approach the full spectral function is decomposed in the con- 
tribution of the QP excitations, and a remainder or background 
part. The latter part is complicated, but does not need to be 
known accurately: it is sufficient to have a functional model 
for the energy-averaged background part to set up a single- 
electron selfconsistency problem that generates the QP exci- 
tations. Such an approach is appealing since it contains the 



well-developed standard Kohn-Sham formulation of DFT as 
a special case, while at the same time emphasis is put on the 
correct description of QPs, in the Landau-Migdal sense [8]. 
Hence, it can provide an improved description of the dynam- 
ics at the Fermi surface. Given the close relation between QP- 
DFT and the Green's function (GF) formulation of many-body 
theory |9|, llfJl . it is natural to employ ab initio calculations 
in the latter formalism to investigate the structure of possible 
QP-DFT functionals. In this respect it is imperative to iden- 
tify which classes of diagrams are responsible for the correct 
description of the QP physics. 

Some previous calculations, based on GF theory, have fo- 
cused on a self-consistent treatment of the self-energy at the 
second order lfTTl[T2i[T3ll for simple atoms and molecules. For 
the atomic binding energies it was found that the bulk of corre- 
lations, beyond Hartree-Fock, are accounted for while signifi- 
cant disagreement with experiment persists for QP properties 
like ionization energies and electron affinities. The formal- 
ism beyond the second-order approximation was taken up in 
Ref. (jj, [lid d H by employing a self-energy of the 
GW type IU9I1 . In this approach, the random phase approx- 
imation (RPA) in the particle-hole (ph) channel is adopted 
to allow for possible collective effects on the atomic excited 
states. The latter are coupled to the sp states by means of dia- 
grams like the last two in Fig.QJc). Two variants of the Go Wo 
formalism were employed in Ref. [ 14] (where the subscript 
"0" indicates that non-dressed propagators are used). In the 
first only the direct terms of the interelectron Coulomb po- 
tential are taken into account. In the second version, also the 
exchange terms are included when diagonalizing the ph space 
[generalized RPA (GRPA)] and in constructing the self-energy 
[generalized GW {GGW)]. Although the exchange terms are 
known to be crucial in order to reproduce the experimentally 
observed Rydberg sequence in the excitation spectrum of neu- 
tral atoms, they were found to worsen the agreement between 
the theoretical and experimental ionization energies fl4ll . 
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In the GW approach the sp states are directly coupled 
with the two-particle-one-hole (2plli) and the two-hole-one- 
particle (2hlp) spaces. However, only partial diagonalizations 
(namely, in the ph subspaces) are performed. This procedure 
unavoidably neglects Pauli correlations with the third particle 
(or hole) outside the subspace. In the case of the GGW ap- 
proach, this leads to a double counting of the second order 
self-energy which must be corrected for explicitly I20I l2~fll . 
We note that simply subtracting the double counted diagram is 
not completely satisfactory here, since it introduces poles with 
negative residues in the self-energy. More important, the inter- 
action between electrons in the two-particle (pp) and two-hole 
(hh) subspaces are neglected altogether in (G)GW. Clearly, it 
is necessary to identify which contributions, beyond GGW, 
are needed to correctly reproduce the QP spectrum. 

In this respect, it is known that highly accurate descrip- 
tions of the QP properties in finite systems can be obtained 
with the algebraic diagrammatic construction (ADC) method 
of Schirmer and co-workers 12211 . The most widely used third- 
order version [ADC(3)] is equivalent to the so-called extended 
2plh Tamm-DancofF (TDA) method 1I23I1 and allows to pre- 
dict ionization energies with an accuracy of 10-20 mH in 
atoms and small molecules. Upon inspection of its diagram- 
matic content, the ADC(3) self-energy is seen to contain all 
diagrams where TDA excitations are exchanged between the 
three propagator lines of the intermediate 2plh or 2hlp prop- 
agation. The TDA excitations are constructed through a diag- 
onalization in either 2p\hox2h\p space, and neglect ground- 
state correlations. However, it is clear that use of TDA leads 
to difficulties for extended systems. In the high-density elec- 
tron gas e.g., the correct plasmon spectrum requires the RPA 
in the ph channel, rather than TDA. 

In order to bridge the gap between the QP description in 
finite and extended systems, it seems therefore necessary to 
develop a formalism where the intermediate excitations in the 
2plh/2hlp propagator are described at the RPA level. This 
can be achieved by a formalism based on employing a set of 
Faddeev equations, as proposed in Ref. [24] and subsequently 
applied to nuclear structure problems 11251 l26l 12711 . In this ap- 
proach the GRPA equations are solved separately in the ph and 
pp/hh subspaces. The resulting polarization and two-particle 
propagators are then coupled through an all-order summation 
that accounts completely for Pauli exchanges in the 2plh/2hlp 
spaces. This Faddeev-RPA (F-RPA) formalism is required if 
one wants to couple propagators at the RPA level or beyond. 
Apart from correctly incorporating Pauli exchange, F-RPA 
takes the explicit inclusion of ground-state correlations into 
account, and can therefore be expected to apply to both finite 
and extended systems. The ADC(3) formalism is recovered 
as an approximation by neglecting ground-state correlations 
in the intermediate excitations (i.e. replacing RPA with TDA 
phonons). 

In this work we consider the Neon atom and apply the F- 
RPA method to a nonrelativistic electronic problem for the 
first time. The relevant features of the F-RPA formalism (also 
extensively treated in Ref. Il24ll ). are introduced in Sect. [TTJ 
The application to the Neon atom is discussed in Sec. [Till 
where we also investigate the separate effects of the ladder 



and ring series on the self-energy, as well as the differences 
between including TDA and RPA phonons. Our findings are 
summarized in Sec. [IV] Some more technical aspects are rele- 
gated to the appendix, where the interested reader can find the 
derivation of the Faddeev expansion for the 2plh/2hlp propa- 
gator, adapted from Ref. Il24ll . In particular, the approach used 
to avoid the multiple-frequency dependence of the Green's 
functions is discussed in App. IA II along with its basic as- 
sumptions. The explicit expressions of the Faddeev kernels 
are given in App. lA 3l Together with Ref. (24]], the appendix 
provides sufficient information for an interested reader to ap- 
ply the formalism. 



H. FORMALISM 

The theoretical framework of the present study is that of 
propagator theory, where the object of interest is the sp prop- 
agator, instead of the many-body wave function. In this paper 
we will employ the convention of summing over repeated in- 
dices, unless specified otherwise. Given a complete orthonor- 
mal basis set of sp states, labeled by a,/3,..., the sp propagator 
can be written in its Lehmann representation as IM UOIl 

gafiiu) = }, — — + >, Z > (1) 

^ 0) - £+ + IT] £f 0)-S k -lTI 

where X" a = <T^ +1 |4|T^> (J/* = <^- 1 |c a |** r » are the 
spectroscopic amplitudes, c a (cl) are the second quantization 
destruction (creation) operators and e+ = E% +1 - E^ (s k = 
E% - ^f- 1 ). In these definitions, PP^ 1 ), PF*" 1 ) are the 
eigenstates, and E% +1 , E^~ l the eigenenergies of the (N ± 1)- 
electron system. Therefore, the poles of the propagator reflect 
the electron affinities and ionization energies. 
The sp propagator solves the Dyson equation 

gap(0j) = g^O}) + g%(uyL* s (uS)g SI3 (Q)) , (2) 

which depends on the irreducible self-energy S*(o>). The lat- 
ter can be written as the sum of two terms 

= ^ajff + J VaA.fjv R/jvA,ySs(^) VySfis > (3) 

where Y, HF represents the Hartree-Fock diagram for the self- 
energy. In Eqs. (0 and (0, g°(u>) is the sp propagator for 
the system of noninteracting electrons, whose Hamiltonian 
contains only the kinetic energy and the electron-nucleus at- 
traction. The V a p t ys represent the antisymmetrized matrix el- 
ements of the interelectron (Coulomb) repulsion. Note that 
in this work we only consider antisymmetrized elements of 
the interaction, hence, our result for the ring summation al- 
ways compare to the generalized GW approach. Equation ([3]) 
introduces the 2p7/!/2/!7p-irreducible propagator R{co), which 
carries the information concerning the coupling of sp states to 
more complex configurations. Both £*(«) and R{to) have a 
perturbative expansion as a power series in the interelectron 
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FIG. 1 : (Color online) a) Diagrammatic expansion of R(u>) in terms 
of the (antisymmetrized) Coulomb interaction and undressed prop- 
agators, b) R(ll>) is related to the self-energy according to Eq. ([3j. 
c) By substituting the diagrams a) in the latter equation, one finds the 
perturbative expansion of the self-energy. 



interaction V. Some of the diagrams appearing in the expan- 
sion of R((o) are depicted in Fig. Q] together with the corre- 
sponding contributions to the self-energy. Note that already at 
zero order in R(ai) (three free lines with no mutual interaction) 
the second order self-energy is generated. 

Different approximations to the self-energy can be con- 
structed by summing particular classes of diagrams. In this 
work we are interested in the summation of rings and ladders, 
through the (G)RPA equations. In order to include such effects 
in R(a>), we first consider the polarization propagator describ- 
ing excited states in the A^-electron system 



H*0 



z 

«*0 



a>-(E»-E») + iTi 
a> + (E»-E»)-iii 



(4) 



and the two-particle propagator, that describes the addi- 
tion/removal of two electrons 



(E»* 2 -E») + in 



y {^\c\clYV N k - 2 ) {^- 2 \ Cf} c \X) 



(5) 



We note that the expansion of R(a>) arises from applying the 
equations of motion to the sp propagator (Q3, which is associ- 
ated to the ground state \¥q ). Hence, all the Green's functions 
appearing in this expansion will also be ground state based, in- 
cluding Eqs. (O and ©. However the latter contain, in their 
Lehmann representations, all the relevant information regard- 
ing the excitation of ph and pp/hh collective modes. The ap- 
proach of Ref. l24ll consists in computing these quantities by 





FIG. 2: (Color online) Diagrammatic equations for the polarization 
(above) and the two-particle (below) propagators in the (G)RPA ap- 
proach. Dashed lines are always antisymmetrized Coulomb matrix 
elements and the full lines represent free (undressed) propagators. 



solving the ring-GRPA and the ladder-RPA equations IK 
which are depicted for propagators in Fig. [2] In the more 
general case of a self-consistent calculation, a fragmented in- 
put propagator can be used and the corresponding dressed 
(G) RPA [D (G)RPA] equations idlH solved [see Eqs. dA2al ) 
and dA2bH . Since the propagators (0 and (O reflect two-body 
correlations, they still have to be coupled to an additional sp 
propagator in order to obtain the corresponding approximation 
for the 2plh and 2hlp components of R(a>). This is achieved 
by solving two separate sets of Faddeev equations. 

Taking the 2plh case as an example, one can split R^ 2plh) (oj) 
in three different components R ( '\a>) (z = 1,2,3) that differ 
from each other by the last pair of lines that interact in their 
diagrammatic expansion, 



1=1,2,3 



arfty,fivA K ' ' 

(6) 

where G (to) is the 2plh propagator for three freely propa- 
gating lines. These components are solutions of the following 
set Faddeev equations B29I1 



R® a ,(a>) = G\^ lllv ,Aco)T w , „,(<u) 

a/3y,(ivA K J apy,fi v A \ / fi'v'A 1 jjf'y"A ' v J 



-(') 



[« 



U) 

fi"v" A" ,fivA 

G° 



(oj) + 

v'M",/iV/t( w ) 



ft v A ,fivA y / 



MP 



v" fi" A" ,fivA 



(«>)] 



(7) 
1,2,3 



where (i,j,k) are cyclic permutations of (1,2,3). The inter- 
action vertices T ( '\cj) contain the couplings of a ph or pp/hh 
collective excitation and a freely propagating line. These are 
given in the Appendix in terms of the polarization (0 and 
two-particle (O propagators. Equations. (O include RPA-like 
phonons and fully describe the resulting energy dependence 
of R(a>). However, they still neglect energy-independent 
contributions-even at low order in the interaction-that also 
correspond to relevant ground-state correlations. The latter 
can be systematically inserted according to 

?(2pl/0 



R i2 R plh \(oj) 



u, 



jjvA,ji'v' A' 



R , ,1 „ U'„ „,„ . 

fi v A ,ft v A v 7 fi v A ,ftvA 



(8) 



where R(oj) is the propagator we employ in Eq. ([3J, R{a>) is 
the one obtained by solving Eqs. (|7), U = I + AU, and I is the 
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FIG. 3: (Color online) Example of one of the diagrams that are 
summed to all orders by means of the Faddeev Eqs. @ (left). The 
corresponding contribution to the self-energy, obtained upon inser- 
tion into Eq. {3}, is also shown (right). 



identity matrix. Following the algebraic diagrammatic con- 
struction method (2~2\ l23ll . the energy independent term AU 
was determined by expanding Eq. ([8]) in terms of the inter- 
action and imposing that it fulfills perturbation theory up to 
first order (corresponding to third order in the self-energy). 
The resulting AU, employed in this work, is the same as in 
Ref. ll23ll and is reported in App. |A 31 for completeness. It 
has been shown that the additional diagrams introduced by 
this correction are required to obtain accurate QP properties. 
Equations. (0 and (O are valid only in the case in which a 
mean-field propagator is used to expand R{ui). This is the case 
of the present work, which employs Hartree-Fock sp propaga- 
tors as input. The derivation of these equations for the gen- 
eral case of a fragmented propagator is given in the appendix. 
More details about the actual implementation of the Faddeev 
formalism to 2plh/2hlp propagation have been presented in 
Ref. 124]]. The calculation of the 2hlp component of R(u>) 
follows completely analogous steps. 

It is important to note that the present formalism includes 
the effects of ph and pp/hh motion to be included simulta- 
neously, while allowing interferences between these modes. 
These excitations are evaluated here at the RPA level and are 
then coupled to each other by solving Eqs. (0. This generates 
diagrams as the one displayed in Fig. [3] with the caveat that 
two phonons are not allowed to propagate at the same time. 
Equations. (0 also assure that Pauli correlations are properly 
taken into account at the 2plh/2hlp level. In addition, one can 
in principle employ dressed sp propagators in these equations 
to generate a self-consistent solution. If we neglect the lad- 
der propagator g"(a>) (0 in this expansion, we are left with 
the ring series alone and the analogous physics ingredients as 
for the generalized GW approach. However, this differs from 
GGW due to the fact that no double counting of the second- 
order self-energy occurs, since the Pauli exchanges between 
the polarization propagator and the third line are properly ac- 
counted for (see Fig. [3]). Alternatively, one can suppress the 
polarization propagator to investigate the effects of pp/hh lad- 
ders alone. 

It is instructive to replace in the above equations all RPA 
phonons with TDA ones; this amounts to allowing only 
forward-propagating diagrams in Fig. [2] and is equivalent to 
separate diagonalisations in the spaces of ph, pp and hh con- 
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TABLE I: Parameters that define the sp basis: radius of the confining 
wall r w (in atomic units) and number of orbits n„ used for different 
partial waves /. The value of c w is always set to 5 a.u.. 



figurations, relative to the HF ground state. It can be shown 
that using these TDA phonons to sum all diagrams of the type 
in Fig.[3]reduces to one single diagonalization in the 2plh or 
2hlp spaces. Therefore, Eqs. and ([8]) with TDA phonons 
lead directly to the "extended" 2plh TDA of Ref. iQ, which 
was later shown to be equivalent to ADC(3) in the general 
ADC framework [22]. The Faddeev expansion formalism of 
Ref. f24ll creates the possibility to go beyond ADC(3) by in- 
cluding RPA phonons. This is more satisfactory in the limit 
of large systems. At the same time, the computational cost 
remains modest since only diagonalizations in the 2p\hj2h\p 
spaces are required. 

Note that complete self-consistency requires the use of 
fragmented (or dressed) propagators in the evaluation of all in- 
gredients leading to the self-energy. This is outside the scope 
of the present paper, but we included partial selfconsistency 
by taking into account the modifications to the HF diagram 
by employing the correlated one-body density matrix and it- 
erating to convergence. This is relatively simple to achieve, 
since the 2plh/2hlp propagator is only evaluated once with 
the input HF propagators. Below we will give results with and 
without this partial selfconsistency at the HF level. 



III. RESULTS 

Calculations have been performed using two different 
model spaces: (1) a standard quantumchemical Gaussian ba- 
sis set, aug-cc-pVTZ for Neon 13011 . with Cartesian repre- 
sentation of the d and / functions; (2) a numerical basis set 
based on HF and subsequent discretization of the continuum, 
to be detailed below. The aug-cc-pVTZ basis set was used 
primarily to check our formalism with the ADC(3) result in 
literature (i.e. HH] , where this basis was employed). The 
HF+continuum basis allows to approach, at least for the ion- 
ization energies, the results for the full sp space (basis set 
limit). 

The HF+continuum is the same discrete model space em- 
ployed previously in Refs. Ill ll 1 1411 . It consists of: (1) Solv- 
ing on a radial grid the HF problem for the neutral atom; (2) 
Adding to this fixed nonlocal HF potential a parabolic poten- 
tial wall of the type U(r) = 6{r — r w )c w (r - r w ) 2 , placed at a 
distance r w of the nucleus. The latter eigenvalue problem has 
a basis of discrete eigenstates. This basis is truncated by spec- 
ifying some largest angular momentum /max and the number 
of virtual states for each value of / < /max- (3) Solve the HF 
problem again, without the potential wall, in this truncated 
discrete space. The resulting basis set is used for the subse- 
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F-TDA F-RPA F-TDAc F-RPAc Expt. 

2p -0.799 -0.791 -0.803 -0.797 (0.94) -0.793 (0.92) 

2s -1.796 -1.787 -1.802 -1.793 (0.90) -1.782(0.85) 

Is -32.126 -32.087 -32.140 -32.102 (0.86) -31.70 

E to t -128.778 -128.772 -128.836 -128.840 -128.928 



F-TDA F-RPA F-TDAc F-RPAc Expt. 

2p -0.807 -0.799 -0.808 -0.801 (0.94) -0.793 (0.92) 

2s -1.802 -1.792 -1.804 -1.795 (0.91) -1.782(0.85) 

Is -32.136 -32.097 -32.142 -32.104 (0.81) -31.70 

Etot -128.863 -128.857 -128.883 -128.888 -128.928 



TABLE II: Results with the aug-cc-pVTZ basis. The first three rows 
list the energies of the main sp fragments below the Fermi level, as 
predicted by different self-energies. F-TDA/F-RPA refers to the Fad- 
deev summation with TDA/RPA phonons, respectively. In all cases 
the self-energy was corrected at third order through Eq. ([8). The 
suffix "c" refers to partial selfconsistency, when the static (HF-type) 
self-energy is consistent with the correlated density matrix. Without 
"c" the pure HF self-energy was taken. In the F-RPAc column the 
strength of the fragment is indicated between brackets. The last row 
is the total electronic binding energy. The experimental values are 
taken from Refs. 133 . 13311 . All energies are in atomic units. 



quent Green's function calculations. 

When a sufficiently large number of states is retained after 
truncation, the final results should approach the basis set limit. 
In particular the results should not depend on the choice of the 
auxiliary confining potential. This was verified in Ref. Il 111 for 
the second-order, and in Ref. flitl for the Go Wo self-energy; in 
these cases the self-energy is sufficiently simple that extensive 
convergence checks can be made for various choices of the 
auxiliary potential. The parameters of the confining wall and 
the number of sp states kept in the basis set was optimized in 
Ref. lITTll . by requiring that the ionization energy is converged 
to about 1 mH for the second-order self-energy. In Ref. IU4ll 
the same choice of basis set was also seen to bring the ioniza- 
tion energy for the Go Wo self-energy near convergence. For 
completeness, the details of this basis are reported in Table J] 
While the self-energy in the present paper is too complicated 
to allow similar convergence checks, it seems safe to assume 
that basis set effects will affect the calculated ionization ener- 
gies by at most 5 mH. 

In Tablellllwe compare, for the aug-cc-pVTZ basis, the ion- 
ization energies of the main single-hole configurations when 
TDA or RPA phonons are employed in the Faddeev construc- 
tion (this is labeled F-TDA and F-RPA, respectively, in the 
table). Note that use of TDA phonons corresponds to the 
usual ADC(3) self-energy. We find that the replacement of 
TDA with RPA phonons provides more screening, leading to 
slightly less bound poles which are shifted towards the exper- 
imental values. This shift increases with binding energy. As 
discussed at the end of Sec. [I] one can include consistency of 
the static part of the self-energy. About eight iterations are 
needed for convergence. This is a nonnegligible correction, 
providing about 5 mH more binding (i.e. larger ionization en- 
ergies) for the valence/subvalence 2p and 2s, 15 mH for the 
deeply bound 1 s, and 60 mH to the total binding energy. Our 
converged result for the Faddeev-TDA self-energy (labeled F- 
TDAc in TableHIli is in good agreement with the ADC(3) value 
for the 2p ionization energy (-0.804 H) quoted in 13111 . as it 
should be. 

The analogous results obtained with the larger 



TABLE III: Results with the HF+continuum basis set from Table U 
See also the caption of Table [TT1 
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FIG. 4: (Color online) Spectral function for the s states in Ne 
obtained with various self-energy approximations. From the top 
down: the second-order (S <2) ) self-energy, the F-RPA(ring), the F- 
RPA(ladder), and the full F-RPA self-energy. The strength is given 
relative to the Hartree-Fock occupation of each shell. Only fragments 
with strength larger than Z > 0.005 are shown. 



HF+continuum basis are given in Table [TTTJ which al- 
lows to assess overall stability and basis set effects. We find 
exactly the same trends as for aug-cc-pVTZ. In particular 
the reduction of ionization energies from the replacement of 
TDA with RPA phonons is almost independent of the basis 
set used, while the effect of including partial consistency 
is roughly halved. Overall, the ionization states are always 
more bound with the larger basis set; while the basis set 
limit could be still more bound than the present results with 
the HF+continuum basis set, it is likely (based on the Go Wo 
extrapolation in Ref. [ 14]) that the difference does not exceed 
5mH. 

As discussed in Sec. |H the the F-RPA self-energy contains 
RPA excitations of both ph type (ring diagrams) and pp/hh 
type (ladder diagrams). It is instructive to analyze their sepa- 
rate contributions to the final ionization energies, in order to 
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Is 2s 2p 

HF -32.77 (1.00) -1.931 (1.00) -0.850(1.00) 

E (2) -31.84 (0.74) -1.736 (0.88) -0.747 (0.91) 

G W Q -31.14 (0.85) -1.774 (0.91) -0.801 (0.94) 

F-RPA(ring) -31.82 (0.73) -1.636 (0.56) -0.730(0.80) 

F-RPA (ladder) -32.04 (0.87) -1.802 (0.95) -0.781 (0.96) 

F-RPA -32.10(0.81) -1.792 (0.91) -0.799(0.94) 

Exp. -31.70 -1.782 (0.85) -0.793 (0.92) 



TABLE IV: Energy (in a.u.) and strength (bracketed numbers) of the 
main fragments in the spectral function of Neon, generated by differ- 
ent self-energies. Results for the HF+continuum basis. Consecutive 
rows refer to: (1) HF; (2) second-order self-energy; (3) Go Wo results 
from Ref. 1 14]; (4) F-RPA self-energy with only ph rings retained; 
(5) F-RPA self-energy with only pp/hh ladders retained; (6) Com- 
plete F-RPA self-energy. In all F-RPA results the self-energy was 
corrected at third order through Eq. l[8}. The static self-energy was 
pure HF (no partial self-consistency). The experimental values are 
taken from Refs. I32L l33fl - 



understand how the F-RPA self-energy is related to the stan- 
dard (G)GW self-energy. Table HVl compares the results for 
the ionization energies, obtained with the second-order self- 
energy, to different approximations for including the ring sum- 
mations. As one can see, the second-order self-energy gener- 
ates an l-\ sp energy of -0.747 mH, which is 46 mH above 
the empirical 2p ionization energy. The Go Wo self-energy, 
which includes the ring summation with only direct Coulomb 
matrix elements, improves this result and brings it close to ex- 
periment. The 2s behaves in a similar way. Unfortunately, 
including the exchange terms of the interelectron repulsion in 
the GG0W3 method turns out to have the opposite effect (the 
2p ionization energy becomes -0.712 H I1 1411 114 111 ) and the 
agreement with experiment is lost. Obviously, GGqWq is too 
simplistic to account for exchange in the ph channel. 

With the F-RPA(ring) self-energy one can go one step fur- 
ther and employ the Faddeev expansion to also force proper 
Pauli exchange correlations in the 2plh/2hlp spaces. As 
shown in Table |IV| this enhances the screening due to the 
exchange interaction terms, leading to even less binding for 
the 2s and 2p. The corrections relative to the second-order 
self-energy can be large (100 mH for the 2s state) and in the 
direction away from the experimental value. We also note that 
the larger shift, in the 2s orbit, is accompanied by an increase 
of the fragmentation (see Fig.|4]and Tab. lIVt . Similar observa- 
tions were also made in Ref. 111411 for other atoms: in general 
ring summations in the direct channel alone bring the quasi- 
hole peaks close to the experiment. This agreement is then 
spoiled as soon as one includes proper exchange terms in the 
self-energy. On the other hand, exchange in the ph channel 
is required to reproduce the correct Rydberg sequence in the 
excitation spectrum of neutral atoms. So further corrections 
must arise from other diagrams, and obviously the summation 
of ladder diagrams can play a relevant role, since these con- 
tribute to the expansion of the self-energy at the same level as 
that of the ring diagrams. 



The result when only including ladder-type RPA phonons 
in the F-RPA self-energy is also shown in Table IIVI One can 
see that pp/hh ladders do actually work in the opposite way 
as the ph channel ring diagrams, and have the same order of 
magnitude with, e.g., a shift of 66 mH for the 2s relative to 
the second-order result. When combined with the ring dia- 
grams in the full F-RPA self-energy, the agreement with ex- 
periment is restored again. Note that the final result cannot be 
obtained by adding the contributions of rings and ladders, but 
depends nontrivially on the interplay between these classes of 
diagrams thereby pointing to significant interference effects. 

With the F-RPA(ring) self-energy, where only the contribu- 
tions of the ph channel are included, the main peaks listed in 
Table [IV] are not only considerably shifted but also strongly 
depleted, e.g. a strength of only 0.56 for the main 2s peak. 
The complete spectral function for the / = strength in Fig. [4] 
shows that the depletion of the main fragment is accompanied 
by strong fragmentation over several states. While correlation 
effects are overestimated in F-RPA(ring), they are suppressed 
in F-RPA(ladder), where only the pp/hh ladders are included 
in the self-energy. In this case one finds a spectral distribu- 
tion closer to the HF one, with a main 2s fragment of strength 
0.95 and less fragmentation than the the second-order self- 
energy. The spectral distribution generated by the complete 
F-RPA self-energy is again a combination of the above ef- 
fects. The strength of the deeply bound Is orbital behaves in 
an analogous way. The strength of the main peak is reduced 
but several satellite levels appear due to the mixing with 2hlp 
configurations. In all the calculations reported in Fig. [4] we 
found a summed I = strength exceeding 0.98 in the interval 
[-40 H, -30 H] which can be associated with the Is orbital, 
and this remains true even in the presence of strong correla- 
tions using the F-RPA(ring) self-energy. Of course, the mix- 
ing with 3h2p configurations, not included in this work, may 
further contribute to the fragmentation pattern in this energy 
region. 



IV. CONLUSIONS AND DISCUSSION 

In conclusion, the electronic self-energy for the Ne atom 
was computed by the F-RPA method which includes - 
simultaneously- the effects of both ring and ladder diagrams. 
This was accomplished by employing an expansion of the 
self-energy based on a set of Faddeev equations. This tech- 
nique was originally proposed for nuclear structure applica- 
tions l24ll and is described in the appendix. At the level of 
the self-energy one sums all diagrams where the three propa- 
gator lines of the intermediate 2plh or 2hlp propagation are 
connected by repeated exchange of RPA excitations in both 
the ph and the pp/hh channel. This differs from the ADC(3) 
formalism in the fact that the exchanged excitations are of 
the RPA type, rather than the TDA type, and therefore take 
ground-state correlations effects into account. The coupling 
to the external points of the self-energy uses the same modi- 
fied vertex as in ADC(3), which must be introduced to include 
consistently all third-order perturbative contributions. 

The resulting main ionization energies in the Neon atom 



are at least of the same quality, and even somewhat improved, 
compared to the ADC(3) result. Note that, numerically, F- 
RPA can be implemented as a diagonalization in Iplh/lhlp 
space implying about the same cost as ADC(3). The present 
study also shows that in localized electronic systems subtle 
cancellations occur between the ring and ladder series. In par- 
ticular, only a combination of the ring and ladder series leads 
to sensible results, as the separate ring series tends to correct 
the second-order result in the wrong direction. 

Since the limit to extended systems requires an RPA treat- 
ment of excitations, the F-RPA method holds promise to 
bridge the gap between an accurate description of quasiparti- 
cles in both finite and extended systems. In particular, the GW 
treatment of the electron gas has been shown to yield excellent 
binding energies, but poor quasiparticle properties 1134 13511 . 
Further progress beyond GW theory requires a consistent in- 
corporation of exchange in the ph channel. The F-RPA tech- 
nique may be highly relevant in this respect. A common 
framework for calculating accurate QP properties in both fi- 
nite and extended systems, is also important for constrain- 
ing functionals in quasiparticle density functional theory (QP- 
DFT) B. 

Finally, complete self-consistency requires sizable compu- 
tational efforts for bases as large as the HF+continuum basis 
used here. It would nevertheless represent an important exten- 
sion of the present work, since it is related to the fulfillment 
of conservation laws ll36ll3~7ll . These issues will be addressed 
in future work. 
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FIG. 5: (Color online) Diagrammatic representation of Eq. {A3)- 
Double lines represent fully dressed sp Green's funcions which, how- 
ever, are restricted to propagate only in one time direction [i.e., only 
one of the two terms on the r.h.s. of Eq. iQJ is retained]. The Faddeev 
Eqs. dA9t and I0 allow for both forward and backward propagation 
of the phonons T ( "\tij) and r (//) (ai) as long as these do not overlap 
in time. For the propagators, time ordereing is asumed with forwad 
propagation in the upward direction. 



of dressing the sp propagator self-consistently. Since the for- 
ward (2plh) and the backward (2hlp) parts of R(co) decouple 
in two analogous sets of equations, it is sufficient to focus on 
the first case alone. 
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1. Multiple frequencies integrals 



APPENDIX A: FADDEEV EXPANSION OF THE 2P1H/2H1P 
PROPAGATOR 



Although only the one-energy (or two-time) part of the 
2plh/2hlp propagator enters the definition of the self energy, 
Eq. ([3]l, a full resummation of all its diagrammatic contribu- 
tions would require to treat explicitly the dependence on three 
separate frequencies, corresponding to the three final lines in 
the expansion of R(co). For example, inserting the RPA ring 
(ladder) series in R(co) implies the propagation of a ph (pp/hh) 
pair of lines both forward and backward in time, while the 
third line remains unaffected. A way out of this situation is 
to solve the Bethe-Salpeter-like equations for the polarization 
and ladder propagators separately and then to couple them to 
the additional line. If it is assumed that different phonons do 
not overlap in time, the three lines in between phonon struc- 
tures will propagate only in one time direction [see figures (|3} 
and ©I- In this situation the integration over several frequen- 
cies can be circumvented following the prescription detailed 
in the next subsection. This approach will be discussed in 
the following for the general case of a fully fragmented prop- 
gator, in order to derive a set of Faddeev equations capable 



We start by considering the effective interactions in the ph 
and pp/hh channels that correspond to Eqs. (|4]i and Q stripped 
of the external legs. In the present work, these are the follow- 
ing two-time objects: 



r (7r 2 Xbi) 



afl,y5> ' 



(Ala) 



ri \ aft) yS V-1 Pa \ SyJ 

VaSfty + / . ; — _ / , : ~ > 

/ — ' co - el + in co + sZ, - w 

V a /3,yS + Vafijty g" %pa -(C0) Vp<r, y 5 (Alb) 

V* V a/3 ) yS V - ! 0)8 V yd ) 
VaSfiy + / . FT — - / . F 

^— < tu - sf+ + irj 4-* co - et~ - 



"I 



where the residues and poles for the ring series are Q^g = 
(^\clc v \^)V^ va and^ = E%-E%. For the ladders, = 
{^ +2 \clcl\X)y,v, a /3 and A-J = W<^-%c v |^>, with 
poles s T n + = E^ +2 -E^ and e£ - = E^-E^ 2 . Equations. dATj 
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solve the ring and ladder RPA equations, respectively 

r aly S (^ = V «WV ( A2a ) 

, n) f du>\ 

r X» = *W ( A2b ) 

i 



To display how the phonons ( lAlat and ( lAlbb enter the ex- 
pansion of R(a>), we perform explicitly the frequency integrals 
for the diagram of Fig. [5] Since it is assumed that the separate 
propagators lines evolve only in one time direction, only the 
forwardgoing {g > (co)) or the backwardgoing (g < (a>)) part of 
Eq. (Q~|i must be included for particles and holes, respectively. 
After some algebra, one obtains 



/d(x>\ da>2 ds d£l 



(A3) 



y, 



X 



fn"" )* fi"" , [tO - 6* - £+, - £+ + £7 - £+ + £" ] fi" ; „ (fif )* 



i v. 



CA + '""V A 

^ari o- 2 ^ 



+■«;/ 

•fllVl 



o- 2 ,/j]Vi 



£0 - (£+ + £+ - £ fc7 ) + JJ7 

1 (<e^ 3 )* <C<*^ fi< , (fit,)* xi\ a^ 2 (a;*')* (KlKl^J *£W* 7 



w - (e 



r- 



L n' 



■ st + ] - ^ + ] [eT 

«2 £3 J L n„ "4 K7 J L ku 



The last term in this expression contains an energy denom- 
inator that involves the simultaneous propagation of two 
phonons. Thus, it will be discarded in accordance with our 
assumptions. It must be stressed that similar terms, with over- 
lapping phonons, imply the explicit contribution of at least 
3p2h/3h2p. A proper treatment of these would require a non 
trivial externsion of the present formalism, which is beyond 
the scope of this paper. 

The remaining part in Eq. iA3i is the relevant contribution 
for our purposes. This has the correct energy dependence of 



a product of denominators that correspond to the intermediate 
steps of propagation. All of these involve configurations that 
have at most 2plh character. Although, ground state corre- 
lations are implicitely included by having already resummed 
the RPA series. Still, this term does not factorize in a prod- 
uct of separate Green's functions due to the summations over 
the fragmentation indices n, and k, [labeling the eigenstates of 
the (N+l)-electron systems]. This is overcome if one defines 
the matrices G 0> {cj), T (1 ' 2 \a>) and T (3) (w), with elements (no 
implicit summation used) 
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-,o> 



a (3npyk y ; iin^vn v Ak,\ 

r (1)> « t n («) 

an a pnpyk y ; [m^vriyAkx x y 



(a>) - 5 na , n 6 n „ nv 5b k 



r (2)> 

Pnpan a yk y \ vn v ^n p Akx 



(CO) = 




Z 



to - (e+ + s^) + it] 



r (3)> „ , (to) 

Q>n a pnpyky \ ftn^vriyAk^ 7 



w + Yj — 



(A + ' n "Y A + *"' 



[&>-<,- < ~ < + ^ ~ < + s kj ] (□";)* I 



Z 



(A4a) 



(A4b) 



(A4c) 



r 



In these definitions, the row and column indices are ordered to 
represent at first two quasiparticle lines and then a quasihole. 
The index T in T®> refer to the line that propagates indepen- 
dently along with the phonon. Using Eqs. (lA4l >. the first term 
on the r.h.s. of Eq. dA3l ) can be written as 



AR (2 ' Ah) ,(<u) = 



V [G 0> ( W )r( 1 >( W )G 0> ( W )r( 3 >( W )G 0> ( W )l 

/ J L J an a pnpyky ; 



(A5) 



/m^vtiyAki 



"a nfi fcy 



Eq. dA51 > generalizes to diagrams involving any number of 
phonon insertions, as long as the terms involving two or more 
simultaneous phonons are dropped. Based on this relation, we 
use the following prescription to avoid performing integrals 
over frequencies. One extends all the Green's functions to ob- 
jects depending not only on the sp basis' indices (a,f3,y) but 
also on the indices labeling quasi-particles and holes («,■ and 
ki). Whether a given argument represents a particle or an hole 
depends on the type of line being propagated. At this point one 
can perform calculations working with only two-time quanti- 
ties. The standard propagator is recovered at the end by sum- 
ming the "extended" one over the quasi-particle/hole indices. 



2. Faddeev expansion 

The 2plh/2hlp propagator that includes the full resumma- 
tion of both the ladder and ring diagrams at the (G)RPA level 
is the solution of the following Bethe-Salpeter-like equation, 



If this equation is solved, a double integration of R{co\ , u>2, C03) 
would yield the two-time propagator R(u>) contributing to 
Eq. (O. However, the numerical solution of Eq. ( IA6t appears 
beyond reach of the present day computers and one needs 
to avoid dealing directly with multiple frequencies integrals. 
The strategy used is to first solve the RPA equations ( lA2al ) 
and ( lA2bb separately. Once this is done it is necessary to re- 
arrange the series (fA6) in such a way that only the resummed 
phonons appear. Following the formalism introduced by Fad- 
deev B29l 13811 . we identify the components R®(co) with the 
three terms between curly brakets in Eq. ( IA6t . By employ- 
ing Eqs. dA2a| ) and flA2b| > one is lead to the following set of 
equations [ 4211 , 



/ds\ 



ds2 dsj, 



x 1 r 

27n "lPm-Mm* 



(to 2 )gy ]y (-coi) (A7) 
(o»i, CO2, coy, s\, S2, S3) 



+R M,»vM> S 2, S3) + W* 1 ' * 2 ' S3)) , i = 1,2, 3 , 



y(k) 



where (i,j,k) are cyclic permutations of (1,2,3) and the inter- 
action vertices F w («i , CO2, 0)3) are given by 



(A8a) 



r afly,iivA (tO\,C02,C03;COi,COs,C0 6 ) = 

= r^ rv/J;1 (w2, W!,w 3 ; w 5 , w 4 ,w 6 ) = 

= 6(coi - co 4 )6(co2 + CO3-CO5- to 6 )g~l(toi)Ijj vA (to 2 + C03) , 



afiy,fivA 

1 



(C0l, C02, CO3; CO4, CO5, CO(,) — 



(A8b) 



R a /5y,HvA(tO\,C02,^3) = (A6) 
= [ga^(C0i)g/s v (c02) ~ gpiiiatdgavfati] gXy{-«0?i) + 

Lp / 3 1 (C02)gy 1 y(-C03)Vfi ll r,y 1 p J ^.Rapa^xitOl, S, C0 2 + CO3 - s) = ^(co 3 - C0 6 )S(cOi + C0 2 - t0 4 - W 5 )^(-W 3 )r^ v (wi + C0 2 ) 

r ds 

+ g aai (0Jl)g yi y(-0J3)V aia - t y lp —R p p a -^ VA {S,C02,C0\+C03-S) 

J lm Finally, we apply the prescription of Sec. lAll and substitute 

1 C ds 1 R{co\,co7, 003) with its extended but two-time version R(co). 

+ -gaa, (coi)g m (co2)V aiJSuP <r J — Rpayj^xKS, toi + to 2 ~ s, m) j . This leads tQ the following set of Faddeev equations which 
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propagate 2plh forward in time, 

an n pn0k y ; ftn^vn v Aki v 7 



^ an a fin.0ky\ [m^vnyAk^^^ ^ an n fin0k y \ fin^vriyAk^^) 



X 



/j'n^v'tt^l'/:^; ftn^vnyAk^ 



(u) + R 



(*) 

pL'nLv'n' v A'k' x \ pm^vriyAkx 



(to) 



i= 1,2,3 . 



(A9) 



Since the full energy dependence is retained in Eq. JA7| >. 
the self-energy corresponding to its solution, R(to\, a>2, <^i), is 
complete up to third order [see Eq. (fj)]. This is no longer the 
case after the reduction to a two-time propagator. In particu- 
lar, the approximation that only forward 2plh propagation is 
allowed between different phonons implies that all diagrams 
with different time propagation of their external lines are ne- 
glected in Eqs. (|A9l >. However, these terms are not energy 
dependent and can be can be reinserted in a systematic way a 
posteriori as in Eq. ©. In the general case, 



^ff/?y,/jvi(w) 



(A10) 



apy;orn a p'nyy ky a'n a p'n^y'k y \ii'n^v'n v A'k^ ' // n^v' n v A' k t ; jivA 



and 

U { Y h) » = 6 a u8B V 6 yl + AI/2f 1A) , (All) 

where the correction AU can be determined by comparison 
with perturbation theory. 



The vertices ( IA4t . that appear in Eqs. ( |A9t , and U^ 2plll) are 
expressed in terms of the fully fragmented propagator. There- 
fore, this approach allows to obtain self-consistent solutions 
of the sp Green's function [25]. Whenever, like in this work, 
only a mean-field propagator is employed as input there ex- 
ist a one-to-one correspondence between the fragmentation 
indices and the sp basis. This is expressed by the relations 
X n a = 8nJ\ - <W) and J/* = 6 k .adaeF, where F represents 
the set of occupied orbits. In this case, it is possible to drop 
one set of indices so that Eqs. (|A9t and ( lAlOb simplify into 
the form © and ©. 



3. Faddeev vertices 

In practical applications, it is worth to note that the poles 
of the free propagator G°(co), Eq. ( IA4ab . do not contribute to 
the kernel of Eqs. ( |A9t . This can be proven by employing the 
closure relations for the RPA problem, in the form obtained 
by extracting the free poles in Eqs. (1A2I >. As an example, for 
the forward poles of the ladder propagator these are 



lim [(w - s* - s*) X (Eq.]A2b]l] => (A 12) 

(X n a '?gyX n ;X? Y^Joj = < + <) = , Vm,n 2 , 

and similarly for other cases. Making use of these relations 
one can derive the following working expression of the ker- 
nels of the 2plh Faddeev equations (no implicit summations 
used) 



J 



[G 0> (w)r (1)> (w)l = \G 0> (u))T (2)> (a))] = (A13a) 

L lan n f3npyk y ; fin^vnyAkx L lf3tipan n yk y ; vtiyfin^Akx 



[G 0> Mr< 3)> ( w )l = (A13b) 

L Jan a pnpyk y ; fm^vriyAkA 

~ K ' kA 2 x<r Kf 1 - [e - - < - £ ^ ][w - (e "" ~ s ^ + iT > ] + if [ C - < - - < - < ] J ' 

After substituting Eq. (lAlOt into ([3), one needs the working expression for the matrix product V U (2plh) (where V is the inter- 
electron interaction). The minimum correction that guaranties to reproduce all third order self-energy diagrams is 

[ VU(2plh) ] Ak = W + ,' r " SA 72 ^ + , (A14) 

L ia; fjn^Ak,, 2[s k + - £+ - £+ ] 

Vojjjiy, {^n^Z) v yi^Siy VaS u vyi {^n^Z) v nMiii 

K + % - < - <] + % - S+ ns ~ <„ ] 

The case of 2/;7p is handled in a completely analogous way along the steps of Sees. (IA It and (IA 2b . After extending 
R(u>\, 0)2, a>j,) to depend on the fragmentation indices (k\Jii,ri), the equivalent of Eq. (|A9t is obtained with the follow- 
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ing definitions of the kernels, 

rtyjx? (jfr)' ( y e^, KT v ggmj^ggjfg | 

E. i < h ^ " % + " " " /7?] + * [ < " % + < ][e ^ " ^ + < ] J ' 

' h '" A 2 Eo- |X-f " ^ " %K» ~ K ~ S V ~ W + if ^ - H a ~ %Mt, ~ % ~ s- k l j ' 

and correction to the external legs, 

17 y"y y"s (y n r y n s\ v 

[ vu(2hlp) ] , ti = w + Vo ™ x * x « ™z ( Ai6) 

Votf!^ (^n^s $ i) v n<Uhv Vafum ^n^e" {^n^j V niton 



It should be pointed out that while the prescription of 
Sec. I A II allows sp lines to propagate only in one time di- 
rection, it allows for backward propagation of the phonons. 
These contributions translate directly into the energy indepen- 
dent terms of Eqs. dA13t and JA15t and are a direct conse- 
quence of the inversion pattern typical of RPA theory. These 
terms have normally a weaker impact than the direct ones on 
the solutions of Eqs. iA9i . However, it is show in Ref. ll24ll 
that they are crucial to guarantee the exact separation of the 
spurious solutions — always introduced by the Faddeev for- 
malism 1 39, 40] — if RPA phonons are used. For the same rea- 



sons, the last terms in curly brackets of Eqs. (IA13b and ( IA15b 
should be dropped whenever Tamm-Dancoff (TDA) phonons 
are propagated. 

The approach followed in this work for solving Eqs. ( IA9I ) 
is to transform them into a matrix representations II24I1 . 
Once this is done, one is left with an eigenvalue prob- 
lem that depends only on the 2plh (2hlp) configurations 
(n, n' ,K) [(k, &',«)]. The spurious states are known ex- 
actly B24I1 and can be projected out analytically to reduce the 
computational load. In any case, they would give vanishing 
contributions to Eq. (0. 
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